† Corresponding author. E-mail:
Project supported by the National Natural Science Foundation of China (Grant No. 11627801) and the Research Foundation of Education Bureau of Hunan Province of China (Grant Nos. 15B083 and 17B090).
The effects of biaxial strain on the electronic structure and thermoelectric properties of monolayer WSe2 have been investigated by using first-principles calculations and the semi-classical Boltzmann transport theory. The electronic band gap decreases under strain, and the band structure near the Fermi level of monolayer WSe2 is modified by the applied biaxial strain. Furthermore, the doping dependence of the thermoelectric properties of n- and p-doped monolayer WSe2 under biaxial strain is estimated. The obtained results show that the power factor of n-doped monolayer WSe2 can be increased by compressive strain while that of p-doping can be increased with tensile strain. Strain engineering thus provides a direct method to control the electronic and thermoelectric properties in these two-dimensional transition metal dichalcogenides materials.
Thermoelectric materials, which can directly and reversibly convert heat into electricity, have potential applications in power generation and refrigeration. The conversion efficiency of thermoelectric materials is described by the figure of merit, ZT, which is defined as ZT = S2σT/κ, where S, σ, and κ are the Seebeck coefficient, electrical conductivity, and thermal conductivity, respectively. A good thermoelectric material must possess a high ZT value. The current researches are focused on enhancing the ZT value by increasing the Seebeck coefficient and electrical conductivity values with reducing the value of thermal conductivity. However, the correlation and coupling among these physical parameters make it extremely difficult to control independently.
The transition-metal dichalcogenides with the formula MX2 (where M = Mo, W; X = S, Se, Te) have found applications in energy storage, sensing, catalysis, and electronic devices for many years.[1–6] These compounds exhibit a “sandwich” type of structure (X–M–X) in which metal atoms (M) are located in between two layers of chalcogen atoms (X).[7] Some crystals of this group, such as MoS2 and WSe2, have been suggested as promising candidates for thermoelectric applications due to the large Seebeck coefficient and low thermal conductivity.[8–14] It was reported that the low-dimensional thermoelectric materials could exhibit much higher ZT values on account of the improved power factor (S2σ) caused by quantum confinement effects.[15,16] Recently, two-dimensional (2D) semiconductor materials formed by transition-metal dichalcogenide layered structures have attracted a great deal of interest due to the fact that they possess enhanced thermoelectric performance compared to the corresponding bulks.[8–14] Compared with other 2D transition metal dichalcogenides, the monolayer WSe2 is found to have an ultralow thermal conductivity, which is one order of magnitude lower than that of MoS2, suggesting that the monolayer WSe2 has great application potential in thermoelectric applications.[17] In previous reports, strain engineering was proved as an effective method to tune the electronic structure and thermoelectric properties of semiconductors.[18–20] Recent studies have shown that the thermoelectric properties of Bi2Te3 and Sb2Te3 can be further improved by the application of strain.[18] It is found that biaxial strain can fine-tune the electronic structure near the Fermi level, and thus it can enhance the transport properties of thermoelectric compounds such as Cu2ZnSnSe4 and SrTiO3.[19,20] Pardo et al. reported that strain can optimize the thermoelectric properties of hole-doped La2NiO4+δ using first-principles calculations as well.[21] It was reported that strain engineering can tune the electronic properties of 2D MX2 (M = Mo, W; X = S, Se, Te),[7–10,22–24] and the strain can be induced easily by growing the transition metal dichalcogenides on flexible substrates. Recently, mechanical strain induced enhancement in the thermoelectric properties have been reported in monolayer ZrS2 and PtSe2 nanosheets as well.[25,26] However, previous theoretical studies mainly focus on the effect of mechanical strains on the electronic properties of the monolayer of MX2 (where M = Mo, W; X = S, Se, Te),[7–10] and few studies have investigated the influence of strain engineering on the thermoelectric properties of WSe2 monolayer. In order to understand the relationship between their electronic structures and thermoelectric properties under strain, explore the effect of strain on their thermoelectric properties, and further enhance the thermoelectric performance of monolayer WSe2, it is necessary to investigate the effect of strain on the electronic structures and thermoelectric properties of monolayer WSe2.
First-principles based on the Vienna ab initio simulation package (VASP)[27] is used to calculate the structural and electronic properties of monolayer WSe2. For the structure optimization and the electronic structure calculations, we employ the Perdew–Burke–Ernzerhof (PBE)[28] generalized gradient approximation (GGA) and projector augment wave (PAW)[29] pseudopotentials, and the cutoff energy of the plane-wave expansion is set at 450 eV. For each monolayer, a vacuum region of 15 Å is added perpendicular to the WSe2 monolayer plane so that the interactions between the periodic images can be neglected. A well-converged Monkhorst–Pack k-point set (13 × 13 × 1) is used for the Brillouin zone sampling. Spin–orbit coupling (SOC) were considered here because it has very important effects on electronic structures and power factor of transition-metal dichalcogenides.[26,30,31]
The transport calculations of the Seebeck coefficient S and electronic conductivity over relaxation time σ/τ are performed through solving the Boltzmann transport equations within the rigid band approximation (RBA)[32] as implemented in the BoltzTraP package.[33] The constant relaxation time approximation,[34] which assumes τ to be energy-independent, is used to estimate the transport properties, and this approximation has successfully predicted the transport properties for many thermoelectric compounds.[18–21,35–38] To get reasonable transport properties, a more dense k-mesh (35 × 35 × 1) was used for WSe2 monolayer in the Monkhorst–Pack scheme, and it can guarantee convergence and obtain accurate carrier group velocities, which is essential for determining the electrical transport properties of WSe2 monolayer.
Bulk WSe2 is crystallized in the hexagonal structure with the P63/mmc space group. The hexagonal structure of WSe2 with the lattice parameters a = b = 3.282 Å and c = 12.96 Å[39] was taken as the starting point for the geometry relaxation. We then built the structure of WSe2 monolayer, starting with the lattice parameters of the bulk relaxed WSe2. The WSe2 monolayer was placed in a large supercell with about 15 Å of vacuum between the periodic replica. The calculated values of a (3.314 Å) are very close to other reported theoretical results.[11] The top and side views of the fully relaxed structural configurations of monolayer WSe2 are shown in Fig.
Band structures of the monolayer WSe2 under different biaxial strains are shown in Fig.
To gain more insight into the changes in the electronic properties of WSe2 monolayer under strains, the partial densities of states (PDOS) of WSe2 monolayer under different applied strains are shown in Fig.
Based on the calculated electronic structure, the Seebeck coefficient S and relaxation time scaled electronic conductivity σ/τ can be obtained using the semi-classical Boltzmann theory in conjunction with rigid band and constant relaxation time approximations. The results for the Seebeck coefficient, electrical conductivity, and power factors of n-doped monolayer WSe2 as a function of the number of electrons per unit cell are shown in Fig.
In order to explore in detail the dependence of thermoelectric properties in p-type doped monolayer WSe2 on strain, the Seebeck coefficient, electrical conductivity, and power factors of p-doped monolayer WSe2 as a function of the number of holes per unit cell are shown in Fig.
In general, 2D transition metal dichalcogenides are n-type semiconductors.[8] The p-type doping of the 2D material can be achieved by chemical doping such as adsorption of small molecules.[42,43] As shown in Figs.
In summary, the influence of biaxial strain on the electronic structure and thermoelectric properties of monolayer WSe2 have been investigated based on first-principles calculations. It is observed that the energy band gap of monolayer WSe2 decreases with the application of biaxial strain, meanwhile, the dispersion of bands near Fermi level of monolayer WSe2 are modified by biaxial strain. It is found that the shift of W d orbital near Fermi level from the projected DOS is caused by the change of W–Se bond length under strain, and it will shrink the band gap and change the shape of the band structure near the Fermi level. The transport properties of n- and p-doped monolayer WSe2 under strain have been estimated based on semi-classical Boltzmann transport theory. Results suggest that compressive strain applied to n-doped WSe2 monolayer is found to be more effective than tensile strain applied onto p-doped WSe2. Meanwhile, it is found that both n-type doping under compressive strain and p-type doping under tensile strain can lead to enhancement of power factor in monolayer WSe2. This research shows that the applied strain is an efficient route for improving the thermoelectric performance of monolayer WSe2 and other 2D transition metal dichalcogenides materials.
[1] | |
[2] | |
[3] | |
[4] | |
[5] | |
[6] | |
[7] | |
[8] | |
[9] | |
[10] | |
[11] | |
[12] | |
[13] | |
[14] | |
[15] | |
[16] | |
[17] | |
[18] | |
[19] | |
[20] | |
[21] | |
[22] | |
[23] | |
[24] | |
[25] | |
[26] | |
[27] | |
[28] | |
[29] | |
[30] | |
[31] | |
[32] | |
[33] | |
[34] | |
[35] | |
[36] | |
[37] | |
[38] | |
[39] | |
[40] | |
[41] | |
[42] | |
[43] |